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ABSTRACT— 

We  are  designing  and  fabricating  arrays  of  subwavelength  SP 
structures  in  metal  films  and  combining  these  structures  with  EO 
polymers  to  understand  the  interaction  of  the  EO  polymer  with 
the  enhanced  electric  fields  of  the  plasmons.  We  have  designed 
and  modeled  resonance  structures  to  maximize  extraordinary 
transmittance.  We  are  making  a  systematic  study  of  SP  structure 
shapes  by  studying  circular  holes,  elliptical  holes  and  slits. 
Initially  we  have  fabricated  and  begun  optical  evaluation  of 
subwavelength  slits  in  gold  films.  Slits  with  32  nm  width  and 
120nm  width  have  been  fabricated  by  a  focused  ion  beam  system. 
Preliminary  optical  transmittance  measurements  in  the  visible  and 
NIR  are  underway  and  will  be  compared  to  theoretical  modeling. 

INTRODUCTION 

Surface  plasmons  (SPs)  are  of  interest  to  a  wide 
spectrum  of  scientists,  ranging  from  physicists,  chemists 
and  materials  scientists  to  biologists  [1],  Renewed  interest 
in  SPs  comes  from  recent  advances  that  allow  metals  to  be 
structured  and  characterized  on  the  nanometer  scale.  This 
in  turn  has  enabled  us  to  control  SP  properties  to  reveal 
new  aspects  of  their  underlying  science  and  to  tailor  them 
for  specific  applications. 

SP-based  photonics  are  being  developed  to  take 
advantage  of  the  larger  information  capacity  of  an  optical 
communications  system  over  an  electronic  one  while 
drastically  reducing  the  physical  size  of  components  over 
other  optical  approaches.  In  an  SP  structure,  the 
electromagnetic  wave,  can  propagate  along  the  metal 
surface  in  the  form  of  an  SP  [2].  As  a  consequence,  the 
fields  can  be  strongly  confined  to  the  metal  surface  with  the 
lateral  dimensions  much  smaller  than  the  wavelength. 
Therefore,  plasmonic  circuits  possess  both  the  capacity  of 
photonics  and  miniaturization  of  electronics,  opening  a  new 
way  for  the  future  applications  [3]. 

SP-based  nanophotonics  or  nanoplasmonics  provide  a 
pathway  to  both  control  and  manipulate  optical  signals  by 
coupling  them  to  coherent  electronic  excitations  near  a 
metal  surface  [1].  This  has  led  to  the  development  of  a 
toolbox  of  various  subwavelength  photonic  components 
such  as  mirrors,  lenses,  plasmonic  crystals,  and  waveguides 
capable  of  manipulating  plasmonic  signals  [4].  Moreover, 


Dipartimento  di  Energetica,  University  of  Rome  La 
Sapienza,  Via  Scarpa  16,  Rome  Italy 
Marco  Centini 

Dipartimento  di  Energetica,  University  of  Rome  La 
Sapienza,  Via  Scarpa  16,  Rome  Italy 

Jason  D.  Fowlkes  and  Ilia  N.  Ivanov 
Center  for  Nanophase  Materials  Sciences 
Oak  Ridge  National  Laboratory,  PO  Box  2008 
Oak  Ridge  TN  37831-6493 

the  plasmonic  properties  of  metallic  nanostructures 
determines  their  conventional  optical  properties  such  as 
reflection,  transmission,  and  absorption.  Plasmonic  crystals 
are  considered  as  the  basis  for  superlenses,  metamaterials 
for  negative  refraction  applications,  as  well  as  nonlinear 
metamaterials  with  enhanced  optical  nonlinearity  [5,  6]. 
Plasmonic  nanostructures  are  becoming  important 
constituents  in  light-emitting  devices  (LEDs,  OLEDs, 
nanolasers),  photodetectors,  and  nanoscale  light  sources  as 
well  as  the  principle  components  for  future  optical  data 
storage,  imaging,  and  sensing.  One  of  the  most  important 
yet  challenging  requirements  for  applications  in  photonics, 
opto-electronics,  and  optical  communications  is  to 
introduce  a  means  to  control  the  optical  properties  of 
plasmonic  devices.  To  directly  influence  surface  plasmons 
and  hence  the  associated  optical  properties,  manipulation  of 
the  dielectric  value  of  the  medium  surrounding  the 
plasmonic  metal  structure  can  be  used  to  change  the 
resonance  condition  for  light-plasmon  coupling.  By 
combining  optically  nonlinear  dielectrics  with  the  metallic 
structure  all-optical  control  can  be  accomplished  by  taking 
advantage  of  the  intense  nature  of  confined  electromagnetic 
radiation  in  nanostmctured  elements  [7,8].  Another 
technique  of  fundamental  importance  in  the  context  of 
controlling  the  SP  resonance  is  the  application  of  an 
external  electric  field  [9].  This  is  especially  significant 
because  both  plasmonic  and  electric  signals  can  be  guided 
in  the  same  metallic  circuitry.  The  applied  voltage  can 
then  be  used  to  exert  direct  control  over  an  electro  optic 
(EO)  dielectric  medium  changing  its  refractive  index  and 
thus  the  SP  modes  at  the  metal/  dielectric  interface. 

This  research  has  two  purposes:  (1)  to  investigate  the 
affect  of  the  enhanced  local  field  of  surface  plasmons  on 
EO  polymers  that  surround  the  plasmonic  structures  as  a 
function  of  hole  shape  and  applied  voltage  and  (2)  develop 
plasmonic  optical  devices  using  plasmonic  structures  and 
EO  materials  for  use  in  photonic  circuits.  Using  an  EO 
polymer  as  a  host  medium  we  will  investigate  how 
nonlinear  optical  materials  and  plasmonics  structures  may 
be  combined  to  create  new  devices.  In  order  to  take 
advantage  of  the  EO  properties  of  the  polymer,  we  are 
designing  and  building  plasmonic  devices  with  resonant 
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properties  in  the  MWIR  and  NIR  wavebands.  We  are 
embedding  the  plasmonics  structures  in  an  EO  polymer  and 
measuring  their  optical  response  in  the  mid  and  near  IR 
waveband.  Our  intent  is  to  make  a  systematic  study  of  hole 
shapes  by  studying  circular  holes,  elliptical  holes  with 
increasing  eccentricity  and  slits.  Because  the  physics  of 
light  transmission  changes  with  this  variation  in  shapes  we 
can  make  a  very  complete  and  systematic  study  of  the 
effects  on  the  EO  polymer  by  the  local  fields  for  this 
plasmonic  structure  [12-16].  For  the  first  phase  of  this 
work  we  have  fabricated  slits  in  gold  films. 

1.  THEORETICAL  MODEL 

The  transmission  through  a  metal  substrate  having  one  or 
more  subwavelength  apertures  is  governed  by  thickness  of 
the  metallic  substrate,  slit  size,  and  inter-slit  distance.  The 
structure  under  consideration  consists  of  gold  layer  whose 
dispersion  profile  is  found  in  Palik’s  handbook  [10],  for 
example.  The  variation  of  these  parameters  is  undertaken  to 
enhance  the  linear  transmittance  and  to  begin  to  appreciate 
the  interaction  between  the  geometrical  and  physical 
mechanisms.  The  calculations  and  results  reported  below 
are  obtained  by  solving  Maxwell’s  equations  using  three 
distinct  approaches:  (1)  Commercially  available  software 
(FEMALB);  (2)  A  home-grown  finite-difference,  time- 
domain  (FDTD)  algorithm,  whose  details  can  be  found  in 
reference  [11];  and  (3)  a  time-domain  fast  Fourier 
transform  beam  propagation  (FFT-BPM)  method  [12,  13, 
14].  Since  the  FDTD  method  is  the  most  widely  adopted 
technique,  we  now  outline  the  FFT-BPM  method  because  it 
is  a  potent  alternative  with  significant  flexibility  [12]. 

In  Gaussian  units,  the  system  of  equations  we  aim  to 
solve  is  as  follows: 

1 0E+4^-ap 
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We  assume  a  right-handed  coordinate  system,  and  p  - 
polarized  (TM)  pump  field  of  the  form: 


The  corresponding  macroscopic  polarization  is  given  by: 


(3). 


The  envelope  functions  E® ,  Ef,  H®,  P‘° ,  PE  contain 
implicit  spatial  dependences  that  for  simplicity  have  been 


omitted.  In  addition,  the  envelope  functions  are  not 
assumed  to  be  slowly  varying,  as  no  approximations  are 
made  when  the  fields  and  polarizations  are  substituted  into 
Maxwell's  equations.  In  other  words,  the  decomposition  of 
the  fields  in  Eqs.(2-3)  is  a  matter  of  mere  convenience. 
Substitution  of  Eqs.(2-3)  into  Eqs.(l)  results  in  a  system  of 
seven  coupled  differential  equations  for  the  fields, 
polarizations,  and  corresponding  currents,  which  in  scaled 
form  are  written  as  follows: 
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We  have  chosen  ~kr  =lpm  as  the  reference 
wavelength,  and  have  adopted  the  following  scaling: 
£  =  z / Ar  and  y  =  y/Ar  are  the  scaled  longitudinal  and 
transverse  coordinates,  respectively;  r  =  ct/  Ar  is  the  time 
in  units  of  the  optical  cycle;  J3  =  2 neb  is  the  scaled  wave 
vector;  (i)  ~  (o  /  (or  is  the  scaled  frequency,  and 
a>r  =  2nd  Ar ,  where  c  is  the  speed  of  light  in  vacuum.  <9, 
is  the  angle  of  incidence  of  the  pump  with  respect  to  the 
normal  direction.  The  linear  dielectric  response  of  silver 
has  been  assumed  to  be  Dmde-like,  as  follows: 

co2 

s((d)  =  1 - ^ —  .  This  function  is  then  fitted  to  the  data 

co 2  +iyco 

at  any  given  wavelength.  For  example,  at  800nm,  the  actual 
data  [Palik]  is  fitted  using  the  set  of  parameters: 
(y 0), co  m)=(0. 06,  6.73) .  The  incident  magnetic  field  was 
assumed  to  be  Gaussian  of  the  form: 

3Cc(5’,<^,r  =  0)  =  Hq  e~^~^  +y  ]/w  ,  with  similar 

expressions  for  the  transverse  and  longitudinal  electric 
fields.  The  equations  are  then  solved  in  the  time  domain 
using  the  well-known,  widely-used  split-step  algorithm.  In 
short,  the  propagation  is  done  in  two  separate  steps:  first  the 
material  interaction  is  taken  into  account,  and  then  the  fields 
are  propagated  in  free  space.  The  equations  are  integrated 
using  a  mix  of  a  predictor-corrector  algorithm  for  the 
material  part,  and  fast  Fourier  transform  for  the  free  space 
propagation  part.  The  main  advantages  that  the  FFT 


approach  offers  over  the  FDTD  method  may  be  summarized 
as  follows:  (1)  calculation  of  the  spatial  derivative  is  much 
more  accurate;  (2)  convergence  is  somewhat  more  rapid,  in 
part  due  to  the  accurate  rendition  of  the  spatial  derivatives, 
thus  allowing  for  larger  temporal  and  spatial  integration 
steps;  (3)  the  method  is  unconditionally  stable,  and  allows 
one  to  choose  a  temporal  integration  step  based  on 
numerical  convergence  rather  than  longitudinal  or 
transverse  discretization  steps. 

A  plot  of  the  calculated  optical  transmittance  as  a  function 
of  wavelength  is  shown  in  Figure  1.  The  blue  plot  is  for 
gold  slits  with  air  in  the  trench.  The  read  plot  is  for  slits 
fdled  with  a  material  with  index  1.67. 


Wavelength  (nm) 


Figure  1.  Transmittance  versus  wavelength  for  120  nm  slit  slits  in  a  gold 
substrate. 

2.  SAMPLE  PREPARATION 

Samples  were  fabricated  by  depositing  gold  films 
supported  on  DSP  Si  wafer,  pyrex  wafer  and  glass  slides 
using  a  Thermionics  e-beam  evaporation  system. 
Thicknesses  were  monitored  by  an  Inficon  quartz  crystal 
thickness  monitor.  Gold  thicknesses  were  either  200  nm  or 
250  nm  depending  on  the  particular  sample  type  being 
fabricated.  For  this  study  we  made  two  sets  of  slits:  (1)  33 
nm  wide  slits  at  a  400  nm  pitch  in  200  nm  gold  films  and 
(2)  120  nm  slits  at  a  720  nm  pitch  in  250  nm  gold  films. 
The  slits  were  written  by  an  FEI  Nova™  600  NanoLab 
DualBeam™  (FIB/SEM)  focused  ion  beam  system  using  a 
Ga  ion  beam.  For  pattern  generation  the  FIB  system  is 
controlled  by  a  Raith  ELPFIY  Quantum  software/hardware 
system.  The  accelerating  voltage  was  30  keV  and  the  ion 
current  was  30  pA.  Single  pass  milling  was  used  in 
conjunction  with  a  dwell  time  of  757  ps  for  the  32  nm  slits 
and  757  ms  for  the  120  nm  slits.  A  16  nm  pixel  size  was 
used  in  both  cases.  The  slits  were  90  pm  long  and  enough 
slits  were  written  in  each  case  to  give  a  written  area  of  90 
pm  x  90  pm.  The  pattern  for  32  nm  wide  slits  is  shown  in 
Figure  2.  The  120  nm  wide  slits  are  shown  in  Figure  3. 


Figure  2.  32nm  line  arrays  (400  nm  pitch)  patterned  in  an  Au  film  (~200 
nm  thick)  on  DSP  Si  wafer.  The  array  size  is  90  pm  x  90  pm. 


Figure  3.  120nm  line  arrays  (700  nm  pitch)  patterned  in  an  Au  film  (~250 
nm  thick)  on  DSP  wafer.  The  array  size  is  90  pm  x  90  pm. 

After  the  samples  were  fabricated  and  the  optical 
transmittance  was  measured  they  were  coated  with  an  EO 
polymer.  The  polymer  is  a  nitrophenylhydrazone  based 
phenoxy  thermoplastic  polymer.  The  material  precursor  is 
mixed  in  a  cyclopentanone  solvent.  For  this  application  is 
solids  concentration  was  21%.  The  solution  is  spin  coated 
to  provide  a  thin  film.  It  was  filtered  through  a  0.2  micron 
filter  just  prior  to  spin  coating.  The  samples  were  spun  at 
1300  RPM  for  20  sec  on  a  standard  photoresist  spin  coater. 
The  sample  was  then  baked  at  100°C  for  10  minutes  and 
then  ramped  to  140°C  over  8  minutes  and  held  at  140°C  for 
25  minutes.  It  was  then  cooled  to  30°C  and  stored  in  N2 
for  24  hrs  before  additional  optical  measurements  were 
made. 

3.  OPTICAL  MEASUREMENTS 

Optical  measurements  were  performed  using  two  Ocean 
Optics  spectrometers  and  a  quartz-halogen  light  source. 
Because  we  are  interested  in  the  optical  response  in  both 
the  NIR  and  the  mid-IR  we  used  both  UV-vis  and  a  NIR 
spectrometer  to  cover  the  waveband  from  400-800  nm  and 
1000-2000  nm.  The  light  source  and  the  spectrometers  use 


fiber  optics  with  collimating  lenses  to  deliver  the  light  to 
the  sample  and  collect  it  from  the  sample.  Figure  4  shows 
the  optical  transmittance  as  a  function  of  wavelength  for 
120  nm  wide  slits  in  a  250  nm  layer  of  Au.  The  plot  is  a 
combination  of  the  results  from  the  UV-vis  spectrometer 
and  the  NIR  spectrometer. 


Figure  4:  Optical  transmittance  as  a  function  of  wavelength  for  120nm 
slits  in  250nm  of  Au 


Figure  5  shows  the  same  sample  after  the  EO  polymer  has 
been  applied.  The  peak  has  been  red  shifted  as  predicted 
by  the  theory. 
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Figure  5:  Optical  transmittance  as  a  function  of  wavelength  for  120nm 
slits  in  250nm  of  Au 
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